A Theoretical Investigation on the Gamma-ray Burst Host 

Galaxies 

Jirong Mao 

Q _ INAF-Osservatorio Astronomico di Brera, Via Bianchi 46, 1-23807, Merate (LC), Italy 

(N' 

K^l Yunnan Observatory, Chinese Academy of Sciences, Kunming, Yunnan Province, 650011, 

^ : China 

Key Laboratory for the Structure and Evolution of Celestial Objects, Chinese Academy of 

Sciences 

\J^ • jirong.mao@brera.inaf.it 

X 



o 






'i 



ABSTRACT 



Long-duration gamma-ray bursts(LGRBs) are believed to be linked with the 
star formation. We adopt a galactic evolution model, in which the star formation 
process inside the virialized dark halo at given redshift can be achieved. In this 
> ■ paper, the gamma-ray burst (GRB) host galaxies are assumed to be the star- 

forming galaxies within the small dark halos. The star formation rates(SFRs) 
OO ■ in the host galaxies of LGRBs at different redshifts have been derived from our 

model with the galactic evolutionary time about a few times of 10^ yr and the 



dark halo mass of about 5 x lO^^Mtn. The related stellar masses, luminosities and 



in 

^ ■ aarK naio mass oi aoout o x iU""iWQ 

Q ■ metallicities of these hosts are estimated as well. We further calculate the X-ray 

and optical absorption of GRB afterglow emission. From our model calculation, 
at higher redshift, the SFR of host galaxy is larger, the absorption in X-ray band 
and optical band of GRB afterglow is stronger, in the condition that the dust and 
metal components are released locally, surrounding the GRB environment. These 
model predictions are compared with the Swift and other observational data. At 
lower redshift 2; < 1, as the merger and interaction of some host galaxies are 
involved, one monolithic physical process is not sufficient to fully explain all 
kinds of observed phenomena. 

Subject headings: dust, extinction — galaxies: evolution — galaxies: star forma- 
tion gamma rays: general 
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Introduction 



Gamma-ray burst (GRB), the most violent explosion cosmic source, has been identified 
as the cosmological event since 1997(van Paradijs et al. 1997; Metzger et al. 1997). Recently 
GRB 090423 has been explored at high redshift above 8(Salvaterra et al. 2009a; Tanvir et 
al. 2009). The long-duration GRB(LGRB) progenitors are proposed to be the massive 
collapsing stars (Woosley 1993; Kumar, Narayan & Johnson 2008). Some long bursts have 
been observed to be associated with supernova events(Hjorth et al. 2003; Stanek et al. 
2003; Malesani et al. 2004; Mazzali et al. 2006; Xu et al. 2008), hence having a common 
star-forming origin(Paczynski 1998). Indeed, long GRBs can be found in the star formation 
galaxies and these galaxies are dominated by the young stellar population(Christensen et 
al. 2004). In general, GRBs favor a metal-poor environment (Fynbo et al. 2006; Kewley 
et al. 2007) and the hosts have low stellar masses (Wiersema et al. 2007). Jakobsson et 
al. (2005) proposed that GRB host galaxies, at least those high redshift (2; > 2) hosts, 
trace the star formation of the universe in an unbiased way. The high global star formation 
rate(SFR) history at redshift larger than 6 (Hopkins & Beacom 2006; Yan et al. 2009) 
indicates the possibility of high-redshift GRB production and the detection of host galaxies. 
From the research of Yiiksel et al. (2008) and Kistler et al. (2009), there could be a link 
from star formation to the GRB production in the high redshift universe, in which the GRB 
luminosity function is involved. Moreover, the evolution of the GRB luminosity function 
has been investigated by Salvaterra et al. (2009b). All of these evidence provide the strong 
clue to study the intrinsic link from SFR to GRB production and the possible evolutionary 
properties of GRBs and their hosts. 

The grains and metals produced by the host galaxy will take effects on the GRB after- 
glow emissions. Thus, the GRB progenitors and their environments can be expressed by the 
absorption features of GRB afterglows. The heavy attenuation in the X-ray band has been 
given in the statistic results from Campana et al. (2010), indicating a dense surrounding 
environment of those GRBs. In the mean while, it is also interesting to understand whether 
this kind of strong attenuation is intrinsically evolved with redshift. On the other hand, the 
characteristics of the corresponding absorption in the optical band are still under debate. 
Although the approximate dust extinction law of GRB host galaxies has been given by Chen, 
Li & Wei (2006) and Li et al. (2008), in order to have an explanation of dust obscuration 
and especially to interpret some X-ray detected but optical faint bursts(so-called dark bursts, 
Akerlof & Swan 2007; Kann et al. 2007; Perley et al. 2009), the physical origin associated 
with the star formation and galactic evolution should be studied in an unified scenario. 

In this paper, we specify one physical model of star-forming and metal-poor galaxies 
being as the hosts of long GRBs, exploiting the physical recipes from Granato et al. (2004). 
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In the general scenario of Granato et al. (2004), at each redshift bin, the SFR and galaxy 
mass in the given dark halo potential well have been calculated, with the effects on the 
kinetic feedback of supernova and central black hole. Under this framework, the different 
evolutionary stages of galaxies and the central black holes with different physical conditions 
have been investigated(e.g., Cirasuolo et al. 2005 about the properties of E/SO galaxies; Lapi 
et al. 2006 about the active galactic nucleus luminosity function; Granato et al. 2006 about 
the submillimeter galaxies). In particular, Mao et al. (2007) calculated the UV luminosities 
and the relative dust attenuation in the star-forming and metal-poor galaxies, Lapi et al 
(2008) estimated the long GRB progenitor rates and redshift distribution. Since the updated 
X-ray/optical observations on the GRB afterglows and host galaxies have been performed 
sequentially by Castro Ceron et al. (2008), Evans et al. (2009), Savaglio, Glazebrook & 
Le Borgne (2009), Levesque et al. (2009a) and Fynbo et al. (2009), in this context, it is 
necessary to further compare some properties calculated by our model with these updated 
observational data. We extend the former calculation from Mao et al. (2007), attempting 
to understand the physical origin of the long GRB production and the GRB environment, 
especially, we reveal that some properties from afterglow emissions and GRB hosts have 
shown the possible intrinsic cosmological evolution. 

Throughout the paper, we adopt cosmological parameters: h = 0.7, Qm = 0.3, and 
^A = 0.7. 



2. Model Predictions 

2.1. Model Review 

In the following we report briefly on some key aspects under the framework of star 
formation in the protogalaxies from our recipes(see also Appendix A of Mao et al. 2007 and 
Lapi et al. 2008 for details). In general, the star formation process and the central black 
hole growth are inside the given virialized dark halo with the mass Mhaio- The cooling gas 
will infall toward the center of the dark halo and form into the stars and galaxy; in the 
mean while it will be heated by the central black hole activity. Thus, the total infalling gas 
Minf = —Mcond — M^^ iucludes two parts: one is the condensation gas toward the center of 
the dark halo Mcond, the other is the gas removed by the central black hole activity M^J . 
The condensation timescale tcond is the maximum between the dynamic timescale and the 
cooling timescale at the halo virial radius. Thus, the cold gas Mcoia = Mcond — (1 — R)M^, — 
M^M — M^^^, where M* = Mcoid/t* is the SFR and R is the fraction of gas transferred to the 
cold component by the evolved stars. Adopting the initial mass function(IMF) by Romano 
et al. (2002), we have R ~ 0.3, M^^ and M^^^ are the feedback from supernova and central 
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black hole respectively. Therefore, with the scaling approximation, we have SFR M^{t) = 
^m/(0)(e"*/*'=°"''— e"'^'''*/*'=°"'')/tco„d(7 — 1/s), where t is the evolutionary time, 7 = 1—R+/3sn, 
PsN is the ratio between star formation feedback by supernova and SFR, s ~ tcond/t* ~ 5. In 
a virialized dark matter halo, the total gas Mj„j(0) is about 18% of the dark halo mass. The 
condensation timescale can be estimated as tcond = 4:XlO^ {{1+ z) / 7)~^'^ {Mhaio/ ^0^^ Mq)^-"^ yr. 
The central black hole quenches the star formation in the halo effectively after the time about 
tBH = 2.5xlO\{l+z)/7)-^-^F{Mhaio/lO^^MQ) yr, where F{x) = 1 for a; > 1 andF(a;) = x-^ 
for X < 1. In other words, the cooling gas inside the virialized dark halo forms into stars and 
galaxy, the star formation process begins and persists with a relatively high rate until a few 
times of 10® yr so that the central seed black hole growth is enough to be a supermassive 
black hole and it shines as a quasar. After that, the central black hole will release the kinetic 
feedback, heat the cold gas, and quench the star formation. 

In this paper, the GRB host is believed to be the young, star-forming galaxies. During 
the forming time about a few times of 10^ yr, the star formation process is violently on- 
going. The masses of these host galaxies are in general less than lO^'^M0(Savaglio et al. 
2009). About these GRB host galaxies, three physical inputs emphasized below can decide 
the whole recipes: (1) redshift z: at different redshifts, the star formation and galactic 
evolution processes are different; (2) dark halo mass M^aio'- as Mao et al. (2007) and Lapi et 
al. (2008) pointed out, the host dark halos in which GRBs occur are relatively small, usually, 
they are less than IQ^'^Mq. In this paper, we select 5 x 10^"*^ M0 as a reference value; this 
is consistent with the simulation results of Courty et al. (2007) and Campisi et al. (2009); 
(3) evolutionary time t: the GRB host galaxy is in the initial stage of the galaxy evolution; 
this initial time t is about a few times of 10^ yr, less than about a few times of 10® yr. This 
value is supported by the observations from Thone et al. (2008) and Han et al. (2010). It 
is noted that at this stage the central black hole seed does not have enough growth to be a 
supermassive black hole; thus, the quasar feedback can be ignored in our calculations. The 
evolutionary time t can be roughly estimated as the starburst time in one starburst galaxy 
as well and the central black hole activity takes negligible effects on the galactic evolution. 
Thus, after these three inputs are given, the SFR in the GRB host galaxies can be decided. 



2.2. Results 

We begin the procedure from the SFR and galactic mass calculation of the hosts, the B- 
band absolute magnitudes of the hosts can be derived from the empirical relation of Savaglio 
et al. (2009). We can also obtain the metallicity distribution by applying the mass-metallicity 
relation of Savaglio et al. (2005). Through the transition from UV band attenuation Auv 
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to the dust absorption A^, we follow the recipes of Mao et al. (2007), in which the results 
of Calzetti et al. (2000) have been adopted, as it is suitable for the high redshift star- 
forming galaxies. We further transfer A^ to X-ray column density Nfj^x using the average 
value obtained by Schady et al. (2010). With the SFR and metallicity properties, the Ay 
distribution can be derived as well. 



2.2.1. Star Formation and Metallicity of GRB Host Galaxies 

As we assume that long GRBs occur inside the young and star-forming galaxies, the 
star formation process plays a key role on the environment of GRB production. With the 
reference values about 5.0 x lO^^M© for a dark halo mass and 5.0 x 10^ yr for the galactic 
evolution time, we have reproduced the SFR at each redshift, the stellar mass of host galaxy 
can be derived as SFR multiplied by the galactic formation timescale. The galactic formation 
timescale can be defined as tg = t ■ ((1 + z) /7)~^-^ {Mhaio/^0^'^ Mq)^-^ yr, which has the same 
index numbers of condensation timescale. We have shown the SFR of GRB host galaxies at 
each redshift in Fig. [1] In this figure, the observational SFR values that spread over a wide 
range, from O.OlMQ/yr to about lOMQ/yr within the redshift bin < 2; < 1, are shown. 
Our model predicts in general that SFR values are larger toward the redshift larger than 1. 

The relation between the SFR and the stellar mass of GRB host galaxy is given in Fig. 
|2l This possible correlation is also mentioned by Savaglio et al. (2009). In our model, we 
illustrate that this correlation at given redshift with a given dark halo mass could be due to 
the growth of the host protogalaxies under the certain galactic formation timescale. However, 
there is no straightforward relation shown by the observational data in Fig. |2J We note that 
within the different dark halos the SFRs and the stellar masses are different. Furthermore, 
some galaxies with relatively larger masses at lower redshift might have experienced twice 
or even more times of starbursts during their lifetimes. Especially, it is found easily in the 
plot that the infrared-selected host galaxies have larger stellar masses. This complicated 
situation indicates that at low redshift, the GRB host galaxies might not have a monolithic 
evolutionary process. The starburst triggered by merging or interaction can happen as well. 
The further discussion about the host galaxies in the low-redshift universe will be given in 
Section 3. 

In the work of Courty et al. (2007), the ratio between SFR and B-band luminosity of 
GRB host galaxy has been investigated with the observational data from Christensen et al. 
(2004). Here, we adopt the observational results by Savaglio et al. (2009). Assuming the 
correlation between SFR and B-band absolute magnitude, we find that the data are well 
described by the scaling relation logSFR = -(0.36 ± 0.01)Mb - (6.72 ± 0.27). Therefore, 
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using this scaling relation we have the B-band absolute magnitude of GRB host galaxies in 
each redshift bin, the results are shown in Fig. [31 It is suggested by Malesani et al. (2009) 
that at higher redshift the GRB host galaxies could be brighter. In our model, we point out 
that this is the result coming from the intrinsic SFR redshift distribution. 

In order to investigate the possible metallicity distribution, in this paper, assuming that 
the mass-metallicity relation and its redshift evolution(Savaglio et al. 2005) are valid for the 
GRB host galaxies as well, after calculation on the stellar masses of GRB hosts, we obtain 
the metallicity evolution as shown in Fig. HI The metallicity values of the host galaxies 
slightly decrease toward the higher redshift. This finding is consistent with that obtained 
by Li (2008). We caution that the metal-poor case is not a necessary input condition in our 
model; thus, the metallicity may not be the essence for GRB production. Further metallicity 
estimations of GRB host galaxies are given in Section 3. 



2.2.2. Afterglow Absorptions 

X-ray Telescope(XRT), one of the instruments onboard the Swift satellite, has supplied 
the important X-ray data in the 0.3-10 KeV band for the GRB research. The Swift-XRT 
analysis has been performed automatically, and the spectral results for Swift-observed GRBs 
have been presented by Evans et al. (2009). Usually the X-ray spectrum is fitted by an 
absorbed power law. Thus, the X-ray photon index and the corresponding neutral hydrogen 
column density Nh,x of each GRB can be achieved. We select each Nh,x value of redshift- 
measured GRB and plot these values in Fig. [51 According to the model described by 
Mao et al. (2007), the UV band absorption Auv = O.35(M,/M0 yr^y-^^{Z/ZQ)^-^ is a 
function of SFR and metallicity. Following the calculation of Mao et al. (2007), we calculate 
from UV attenuation to E{B - V) using E{B - V) = Auv/H by Calzetti et al. (2000). 
The results are in agreement with the observations (see Fig. 2 of Mao et al. 2007). With 
R^ = 3.1, we obtained the dust attenuation A^. With the data observed by Swift-XRT 
and Swift-UV/Optical Telescope (UVOT), Schady et al. (2007 & 2010) modeled the spectral 
energy distributions and derived the ratio between Nh,x and Ay. As the ratios derived from 
Schady et al. (2010) might be varied with the redshift, assuming the linear relation between 
redshift and the ratio with a logarithmic scale, we perform the linear regression on the data 
and obtain the optimized relation as log{NH,x/Ay 10'^^cm~'^) = 1.24:log{l + z) + 0.79 with 
the average standard deviation 0.37. Then we use this relation to transfer Ay to Nh^x and 
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compare the results to the X-ray absorption datalll in Fig. [5] panel (a). 

However, we note that the selection effects are included in the results of Evans et 
al.(2009) and Campana et al.(2010). As pointed out by Campana et al. (2010), at high 
redshift larger than 4, the intrinsic X-ray emissions suffer lower absorption; thus, the X-ray 
afterglows with low X-ray column densities are hard to be identified by Swift-XRT. On the 
other hand, as we use the data from Schady et al. (2010), although it was claimed that 
generally the selection effects on the distribution of host column densities are not significant, 
in our paper, in order to investigate the possible selection effects on our results, first, we 
check the possible Nh,x — z relation and the A^ — z relation respectively from the data 
of Schady et al. (2010). We see that the correlation between Nh,x and redshift has the 
efficient r = 0.67 with null hypotheses 0.0006, this relation could be due to the selection 
effect mentioned above, but we do not find any possible relation between A^ and redshift. 
As Nh,x increases with redshift, Nh,x/A^ also increases with redshift. Second, aiming to 
avoid this selection effect, we use the average value < Nu^xjA-u >= 3.3 x lO^^cm"^ given by 
Schady et al. (2010) to transfer A^ to Nh,x again. We plot the results in Fig. [5] panel (b). 
By using the mean value of N^^x/Ay, the selection effect can be effectively depressed. 

After the depression of selection effect, our model results still show a slight trend of 
X-ray absorption evolution. This evolution trend may be intrinsic. From our model, we 
see that the X-ray absorption is originally from the SFR. SFR has the redshift evolution as 
SFR ~ (1 -|- 2;)^''^^. Under the assumption of solar metallicity, we have the intrinsic X-ray 
attenuation Nh,x ~ (1 + z)^'^^. Therefore, we conclude that the SFR redshift evolution is 
the dominant reason for the X-ray attenuation evolution shown in Fig. [5] (b). If we use 
the linear relation between Nh,x/A^ and redshift, meaning the possible selection effects are 
included, we have the final results shown in Fig [5] (a). We see that the intrinsic evolution 
plus the selection effects can fit the observational data of Evans et al. (2009) and Campana 
et al. (2010) well. 

Through the analysis above we clearly see, that the final results of GRB X-ray absorption 
are the calculations of intrinsic SFR redshift evolution, modified by the variation between 
Nh,x/Av and redshift. The later could be due to the selection effect. From Fig. El we see 
that the observational data have large scatter. On the other hand, our model provides the 
different values under the different dark halo masses and evolutionary time. Therefore, we 
also conclude that the large absorption is due to the longer galactic evolution time within 



^As the XRT spectra are processed by the standard software XSPEC, in which the metalhcity is fixed 
as the solar value, thus the observational data and the calculations of Nh,x are all uniformed by the solar 
metallicity. 



the massive dark halo, while the small attenuation is due to the shorter galactic evolution 
time within the smaller dark halo. From the theoretical point of view, we confirm that the 
absorption is from the local environment of GRB, as suggested by Campana et al. (2010), 
Nardini et al. (2009), and Zheng et al. (2009) from the data analysis. 

The quantities of neutral gas in the host galaxies can be obtained from the optical 
spectra. By the measurement of Lya absorption, Fynbo et al. (2009) have established one 
sample in which 33 values of neutral hydrogen column density NH^opt are derived. The range 
of these values is from 10^'^ cm~^ to 10^^cm~^(see Fig. 10 of Fynbo et al. 2009), while the 
true distribution of Nn^opt niay extend to the higher column densities. On the other hand, 
the damped Lya system with the neutral hydrogen number exceeding 2 x 10^°cm~^ has the 
possibility of star formation to form a protogalaxy (see the simulations by Pontzen et al. 
2009 recently); thus, it could be the GRB host. But the thin cloud with the smaller neutral 
hydrogen values might be intervening along the line of sight between the observer and the 
GRB place; this kind of thin cloud with the column density less than ~ 10^°cm~^ might 
not be related to the GRB host. Also, in this paper, we assume that GRB hosts are rich in 
neutral gas. Therefore, we only select the Nfj^opt values larger than 10^^cm~^ and compare 
them with the corresponding X-ray absorption Nh,x values. We find the relation between 
X-ray absorption and optical neutral gas shown in Fig. |6l logNH,x = {0.4:9±0.04:)logNH,opt + 
(11.3 ± 0.9). The linear correlation coefficient is 0.58 with the probability 0.001. Through 
this weak relationship, it is likely to find the trace of possible cosmic evolution of neutral 
gas Nn^opt, similar to the evolution of X-ray Nh^^ in Fig. O Suppose that GRBs are the 
unbiased tracers of star formation at high redshift, this possible Nfj^opt distribution may 
give an interpretation to the observations of HI gas evolution by Prochaska & Wolfe (2009). 
However, we caution that the relation between X-ray absorption and optical neutral gas may 
have larger uncertainties, due to the limited redshift range from 2 to 3. A complete sample 
is required to investigate this relation in the future. 

From our model, we see that dust absorption Ay is the function of SFR and metallicity. 
Combining the effects of both SFR and metallicity, we obtain the redshift distribution of dust 
absorption shown in Fig. [71 We see that the Av values from our model are slightly increasing 
with redshift. From the data of Schady et al. (2010), we do not find any prominent evidence 
of Ay variation. But it was claimed by Kann et al. (2007) that the Ay value decreases with 
increasing redshift. Here, after the comparison between the two data sets given by Schady 
et al. (2010) and Kann et al. (2007), we find that the two data sets are not consistent with 
each other: some Ay values of the same burst have large difference. The details are listed in 
the caption to Fig. [71 At high redshift, the values of Kann et al. (2007) are lower than those 
of Schady et al. (2010), while at low redshift, the values of Kann et al. (2007) are larger 
than those of Schady et al. (2010). In general, low-mass stars take a long time to evolve 
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to asymptotic giant branch (AGB) phase and to produce the dust; thus, AGB population 
only dominate the dust production at local universe. It is suggested that at high-redshift 
the dust factory is supernova explosion. However, recently, AGB population has been found 
to be the source of dust production at high redshift universe (Valiante et al. 2009): dust 
production is mainly from supernova at the beginning of the evolution, but for the time 
larger than 3x10^ yr, the dust contribution from AGB stars increases. From the calculation 
of Valiante et al. (2009), we see that at time 1.0 x 10^ yr, AGB dust production is still 10 
times lower than that of supernova. Thus, in our paper, including the AGB dust production, 
the dust extinction A^ is added by a factor of 8%. While at the time 3 x 10^ yr, the AGB 
dust production is as same as supernova dust production. Thus, we estimate from our model 
that the total dust extinction A^ is 1.7 times of the original value in which only supernova 
production is included. Therefore, we clearly see that both AGB and supernova are the 
origin of dust production at the late evolution time larger than 10^ yr. Finally, we cannot 
ignore the selection effect: at high redshift, the GRBs and their hosts with high absorption 
values are difficult to be detected by optical telescopes. 



3. Discussions 

Under the framework of galaxy formation scenario, Lapi et al. (2008) predicted the GRB 
progenitor rate and redshift distribution. In this paper, without the information of GRB 
rates and the cosmological star formation density, we attempt to reveal some properties of 
GRBs and their host galaxies, which have intrinsic redshift distributions. The distributions 
of these properties with redshift are found to be originally from the star formation in the 
star-forming galaxies. Given a proper galactic evolutionary time and a reasonable dark halo 
mass, the final results can be obtained by the model calculation. These results are compared 
with all kinds of observational data. 

At high redshift, the GRB host galaxy has a plenty of neutral gas, suffering much 
violent star formation. After the short-time stellar evolution phase, the metal and dust are 
released by massive stars; thus, the optical and X-ray GRB emissions will have a strong 
attenuation locally at high redshift. The star formation activity, evolving from relatively 
massive hosts at high redshift to dwarf galaxies at low redshift, is a coincidence with the 
so-called downsizing scenario(e.g.. Heavens et al. 2004). However, at lower redshift, the 
situation turns to be more complicated. From the morphological statistics by Conselice et 
al. (2005) and Wainwright et al. (2007), the GRB hosts present a broad diversity of galaxy 
types. About 1/3 host galaxies in the sample of Savaglio et al. (2009) are mergers, while 
in our model the merging and interaction processes are not taken into account. In fact. 
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Conselice et al. (2005) found that the GRB hosts at z > 1 are different from those aX z < 1 
in terms of hght concentration and the morphological size. Through the study of galaxy mass 
distribution, GRB hosts tracing star formation might be biased at low redshift(Kocevski et 
al. 2009). It is also complex that the hosts at z < 1 are not representative of the general 
galaxy population(Levesque et al. 2009a). Thus, the properties of these low-redshift GRB 
hosts presented in this paper could not be reproduced by any monolithic process. At least, 
some low-redshift galaxies may undergo multiple star-forming processes during their whole 
lifetimes. GRB production can be accompanied with any single starburst event. 

From the analysis in this paper, we see that the absorption of GRB X-ray and optical 
emissions is relatively strong. The strong intrinsic attenuation of GRB host galaxies may 
produce some dark bursts, defined by the index Pox < 0.5, where Pox is the flux density 
ratio between optical and X-ray bands(Jakobsson et al. 2004). Rol et al. (2005) proposed 
several extinction origins from their preliminary results. From our calculations, we see 
that the heavy attenuation may occur due to the following three possibilities: (1) the local 
environment of the host is metal-enriched, metallicity is higher, and/or, the host galaxy in 
the massive dark halo larger than IO^^Mq may have strong absorption. For example, at 
redshift 2.5, Z = I.OZq, halo mass M^aio = 5.0 x lO^^M©, after the galactic evolving time 
1.0 X 10^ yr, we have dust extinction A^ = 1.0 and the corresponding X-ray absorption 
Nh,x = 4.7 X lO^^cm"^; (2) the dust and metals surrounding the GRB in the host galaxy are 
distributed in an inhomogeneous way; there could be heavy absorption through the line of 
sight, but in other directions the absorption is slight. Also, in our model, we assume that the 
Ay and Nh^x are measured locally and do not change significantly if the dust and gas extend 
out to a few tens to hundreds of pc from the burst(Perna & Lazzati 2002, D'Elia et al. 2009); 
however, suppose the observed optical extinction is due to the grain absorption far beyond 
this local region of GRB, the A^-Nh^x correlation obtained by Schady et al. (2007, 2010) 
may be invalid and our calculations are strongly biased; (3) as mentioned in Section 2.2.2, 
the dust produced by the AGB population at high redshift should be taken into account. 

In order to further understand the metal production of GRB environment, we roughly 
re-estimate the metallicity of GRB hosts under our framework. The mass of metal M^etai = 
SFR ■ f ■ fdep ■ Mdust/Mstar, where f is the ratio of massive stars to all stars, and fdep is the 
ratio of mental converted from the dust. / = 0.47 is the case for the stars with the mass 
larger than 2Mq by our adopted IMF, fdep = 1.0 means that all the dust can be transferred 
to metals. Metallicity is defined hj Z = Mmetai/Mgas- From the SFR calculated by Granato 
et al. (2004) and Mao et al. (2007), as an example, at redshift 6, we obtain the metallicity 
as Z ^ 2.75 X 10^^ (Mdust/Mstar)- If we take a supernova with the dust production of 10~^ 
solar mass(Pozzo et al. 2004), we obtain the upper limit of metallicity Z ~ 10~^ Zq, which 
is lower than the measurement(Z > O.O2Z0) of GRB 050904(Campana et al. 2007). If we 
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take the dust mass 0.08-0.3 solar mass per primordial massive supernova(Todini & Ferrara 
2001), we have the result which is consistent with the observation. The estimation values of 
Population I/II metallicity are lower than the observational values at high redshift, meaning 
that the imprints from some primordial objects(Kawai et al. 2006), such as pop III stars 
and mini-quasars, have to be included in the possible cosmic evolution properties of these 
GRB host galaxies. According to this estimation, the metal-enriched environment of GRB 
host galaxy naturally gives the reason of the strong attenuation in X-ray and optical band 
measurements. In our model, the initial galactic evolutionary time of about 10^ yr of host 
galaxies is given, but the corresponding metallicity about Z ~ 0.3ZQ(Lapi et al. 2008) is not 
a necessary condition, as mentioned by Levesque et al. (2009b) that low metallicity may not 
be required for a relativistic explosion. With our model, the massive dark halo above lO^^M© 
can host the GRB galaxy in which the metallicity is relatively high, although most GRB 
host galaxies are inside the dark halos with the masses less than IO'^^Mq. On the other hand, 
a host galaxy with a top-heavy IMF, meaning that much more massive stars are involved, 
can produce more metals in relatively short time during the galactic evolution phase. For 
instance, the Wolf-Rayet star with a mass of 8OM0 and initial metallicity Z = 0.001 has the 
possibility to self-enrich the HII region (Kroger, Hensler & Freyer 2006) and to produce the 
GRB event (Eldridge et al. 2006). 

In this paper, we have calculated the SFR, galactic mass, and metallicity of the GRB 
host galaxies. The absorption variations with redshift in the X-ray and optical bands are pre- 
sented as well. Some selection effects have been taken into account through our calculation. 
Other observational biases should also be considered. All the redshift measurements come 
from the optical observations so that some optical-faint GRBs and host galaxies are ignored. 
Moreover, at high redshift only most luminous galaxies with high SFRs can be detected, 
indicating that some low-luminosity cases are not included. However, our calculations come 
from the intrinsic star formation of GRB host galaxies. Thus, due to all these selection effects 
and observational bias mentioned in the paper, the intrinsic properties of GRB afterglows 
and the hosts by the model calculations have some differences to those from observations. 
As SFR evolution plays a dominant role in the calculations, compared to the situation at 
low redshift, in general, star formation in the metal-poor environment at high redshift may 
provide more powerful GRB explosion. Therefore, although the effective threshold is given 
by Kistler et al. (2009), we speculate that improving the sensitivity of detectors on the 
high-energy telescopes is not strongly useful to catch more high-redshift but faint GRBs, 
since low-energy-released GRBs are almost absent in the high-redshift universe. 
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Fig. 1. — SFRs of GRB host galaxies. The dash-dotted hne, the sohd hne and the dashed 
hne denote the prediction results by the model under the condition of dark halo mass 5.0 x 
IO^^Mq, their corresponding galactic evolutionary times are 2.5 x 10^ yr, 5.0 x 10^ yr and 
7.5 X 10^ yr respectively. The dotted hne and the dash-double-dotted line are the lower and 
upper limits given from the model, with the two set parameters (halo mass 1.0 x IO^^Mq, 
galactic timescale 1.0 x 10^ yr and halo mass 5.0 x IO^^Mq, galactic timescale 1.0 x 10^ yr 
respectively). The observational SFR values taken from Savaglio et al. (2009) are shown as 
dots, from Castro Ceron et al. (2008) as squares, and from Levesque et al. (2009a) as stars. 
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Fig. 2. — Stellar mass as a function of SFR. The model prediction range is dependent on the 
adopted galactic formation timescale. The model predicted solid and dash lines are in the 
condition of redshift 1.0, with the halo masses of 5 x IO^^Mq and 1 x IO^^Mq, respectively. 
The galactic timescale ranges are tg ~ ]^q7.o-8.5 y^ ^^j, ^^^e results denoted by the solid 
line, and tg ^ iQ7.7-9.0y^ ^^j, ^j-^g results denoted by the dashed line. The prediction at 
redshift 0.1 is denoted as the dash-dotted line, the result is under the condition of halo 
mass 1.0 X IO^^Mq and the galactic timescale range tg ~ iq7-5-9.2 ^^_ ^^ redshift 3.0, we 
get the results denoted by the dotted line, under the condition of 5.0 x IO^^Mq and the 
timescale range tg ~ xo6-7-8.2 y^^ r^^^e near-infrared-selected GRB hosts(denoted as squares; 
Castro Ceron et al. 2008) have larger mass values than the hosts selected by the optical 
observations(denoted as dots; Savaglio et al. 2009). We specify the data at redshift larger 
than 2, shown as stars in the plot. 
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Fig. 3. — B band absolute magnitudes of GRB host galaxies. All the model predictions are 
the same as those in Fig. 1. The observational data from Savaglio et al. (2009) are marked 
by dots while those from Levesque et al. (2009a) are marked by squares. 
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Fig. 4. — Metallicity distribution with the normahzation 12 + log{0/H) = IoqIZ/Zq) +8.69. 
The model predictions are the same as those in Fig. 1. The average values from Savaglio et 
al. (2009) are denoted by dots, and the data from Levesque et al. (2009a) are marked by 
triangles. The metallicity values of GRB-damped Lya systems (Savaglio 2006) are presented 
by squares. 
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Fig. 5. — X-ray absorption Nh,x distribution with redshift. All the model predictions are 
the same as those in Fig. 1. The observational data are taken from Evans et al. (2009) and 
Campana et al. (2010). In the plot, the dark bursts identified by Perley et al. (2009) and 
Zheng et al. (2009) are labeled as stars. The results from the model shown in the upper 
panel (a) is the Nh^x distribution with the possible selection effects. The model results with 

thp Hfinrpssinn nf sfilpntinn pffpnts are shown in thp Inwpr nanpl fbV 
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Fig. 6. — Correlation between X-ray absorption Nh,x and neutral hydrogen column density 
NH,opt- These bursts (redshifts) are as follows: GRB 050319 (3.240), GRB 050401 (2.899), 
GRB 050730 (3.968), GRB 050820A (2.612), GRB 050922C (2.198), GRB 060115 (3.533), 
GRB 060206 (4.048), GRB 060210 (3.913), GRB 060707 (3.425), GRB 060714 (2.711), GRB 
060906 (3.686), GRB 060926 (3.206), GRB 060927 (5.464), GRB 061110B (3.433), GRB 
070110 (2.351), GRB 070506 (2.308), GRB 070611 (2.041), GRB 070721B (3.628), GRB 
070802 (2.455), GRB 071031 (2.692), GRB 080210 (2.641), GRB 080413A (2.433), GRB 
080603B (2.690), GRB 080607 (3.037), GRB 080721 (2.591) and GRB 080804 (2.205). 
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Fig. 7. — Dust absorption A^ distribution with redshift. All the model predictions are 
the same as those in Fig. 1. The observational data taken from Kann et al. (2006) are 
denoted as squares, from Kann et al. (2007) as '*', and from Schady et al. (2010) as dots. 
We note some different A^j values of the same burst measured by Kann et al. (2007) and 
Schady et al. (2010): GRB 050730(z=3.967) A^l = 0.10 ± 0.02 and A^ = 0.23 ± 0.02; 
GRB 060206(z=4.048) A^l = 0.01 ± 0.02 and A^2 < 0.18; GRB 060526(z=3.211) A^l = 
0.02±0.03 and A^2 < 0.16; GRB 050820A(z=2.612) A^l = 0.07±0.01 and A^2 = 0.14±0.03; 
GRB 050922C(z=2.198) A^l = 0.01 ± 0.01 and A^2 = 0.07 ±0.02; GRB 050525A(z=0.606) 
A^l = 0.32 ± 0.20 and A^2 = 0.16 ± 0.02; GRB 060904B(z=0.703) A^l = 0.44 ± 0.05 and 
Ay2 = 0.12 ± 0.05. A^l and A^2 are the values given by Kann et al. (2007) and Schady et 
al. (2010) respectively. In general, the values measured by Kann et al. (2007) have larger 
error bars than those measured by Schady et al.(2010). 



